Search for: All records

Abstract Skin coloration and patterning play a key role in animal survival and reproduction. As a result, color phenotypes have generated intense research interest. In aposematic species, color phenotypes can be important in avoiding predation and in mate choice. However, we still know little about the underlying genetic mechanisms of color production, particularly outside of a few model organisms. Here we seek to understand the genetic mechanisms underlying the production of different colors and how these undergo shifting expression patterns throughout development. To answer this, we examine gene expression of two different color patches(yellow and green) in a developmental time series from young tadpoles through adults in the poison frogOophaga pumilio.We identified six genes that were differentially expressed between color patches in every developmental stage (casq1, hand2, myh8, prva, tbx3,andzic1).Of these,hand2, myh8, tbx3,andzic1have either been identified or implicated as important in coloration in other taxa.Casq1andprvabuffer Ca²⁺and are a Ca²⁺transporter, respectively, and may play a role in preventing autotoxicity to pumiliotoxins, which inhibit Ca²⁺-ATPase activity. We identify further candidate genes (e.g.,adh, aldh1a2, asip, lef1, mc1r, tyr, tyrp1, xdh), and identify a suite of hub genes that likely play a key role in integumental reorganization during development (e.g., collagen type I–IV genes, lysyl oxidases) which may also affect coloration via structural organization of chromatophores that contribute to color and pattern. Overall, we identify the putative role of a suite of candidate genes in the production of different color types in a polytypic, aposematic species. 

ABSTRACT The skin microbiome is an essential line of host defense against pathogens, yet our understanding of microbial communities and how they change when hosts become infected is limited. We investigated skin microbial composition in three North American bat species (Myotis lucifugus,Eptesicus fuscus, andPerimyotis subflavus) that have been impacted by the infectious disease, white-nose syndrome, caused by an invasive fungal pathogen,Pseudogymnoascus destructans. We compared bacterial and fungal composition from 154 skin swab samples and 70 environmental samples using a targeted 16S rRNA and internal transcribed spacer amplicon approach. We found that forM. lucifugus, a species that experiences high mortality from white-nose syndrome, bacterial microbiome diversity was dramatically lower whenP. destructanswas present. Key bacterial families—including those potentially involved in pathogen defense—significantly differed in abundance in bats infected withP. destructanscompared to uninfected bats. However, skin bacterial diversity was not lower inE. fuscusorP. subflavuswhenP. destructanswas present despite populations of the latter species declining sharply from white-nose syndrome. The fungal species present on bats substantially overlapped with the fungal taxa present in the environment at the site where the bat was sampled, but fungal community composition was unaffected by the presence ofP. destructansfor any of the three bat species. This species-specific alteration in bat skin bacterial microbiomes after pathogen invasion may suggest a mechanism for the severity of white-nose syndrome inM. lucifugusbut not for other bat species impacted by the disease. IMPORTANCEInherent complexities in the composition of microbiomes can often preclude investigations of microbe-associated diseases. Instead of single organisms being associated with disease, community characteristics may be more relevant. Longitudinal microbiome studies of the same individual bats as pathogens arrive and infect a population are the ideal experiment but remain logistically challenging; therefore, investigations like our approach that are able to correlate invasive pathogens to alterations within a microbiome may be the next best alternative. The results of this study potentially suggest that microbiome-host interactions may determine the likelihood of infection. However, the contrasting relationship between Pd and the bacterial microbiomes ofMyotis lucifugusandPerimyotis subflavusindicate that we are just beginning to understand how the bat microbiome interacts with a fungal invader such as Pd. 

Abstract Background Phylogenomic approaches have great power to reconstruct evolutionary histories, however they rely on multi-step processes in which each stage has the potential to affect the accuracy of the final result. Many studies have empirically tested and established methodology for resolving robust phylogenies, including selecting appropriate evolutionary models, identifying orthologs, or isolating partitions with strong phylogenetic signal. However, few have investigated errors that may be initiated at earlier stages of the analysis. Biases introduced during the generation of the phylogenomic dataset itself could produce downstream effects on analyses of evolutionary history. Transcriptomes are widely used in phylogenomics studies, though there is little understanding of how a poor-quality assembly of these datasets could impact the accuracy of phylogenomic hypotheses. Here we examined how transcriptome assembly quality affects phylogenomic inferences by creating independent datasets from the same input data representing high-quality and low-quality transcriptome assembly outcomes. Results By studying the performance of phylogenomic datasets derived from alternative high- and low-quality assembly inputs in a controlled experiment, we show that high-quality transcriptomes produce richer phylogenomic datasets with a greater number of unique partitions than low-quality assemblies. High-quality assemblies also give rise to partitions that have lower alignment ambiguity and less compositional bias. In addition, high-quality partitions hold stronger phylogenetic signal than their low-quality transcriptome assembly counterparts in both concatenation- and coalescent-based analyses. Conclusions Our findings demonstrate the importance of transcriptome assembly quality in phylogenomic analyses and suggest that a portion of the uncertainty observed in such studies could be alleviated at the assembly stage. 

Abstract Iridescence is widespread in the living world, occurring in organisms as diverse as bacteria, plants, and animals. Yet, compared to pigment-based forms of coloration, we know surprisingly little about the developmental and molecular bases of the structural colors that give rise to iridescence. Birds display a rich diversity of iridescent structural colors that are produced in feathers by the arrangement of melanin-containing organelles called melanosomes into nanoscale configurations, but how these often unusually shaped melanosomes form, or how they are arranged into highly organized nanostructures, remains largely unknown. Here, we use functional genomics to explore the developmental basis of iridescent plumage using superb starlings (Lamprotornis superbus), which produce both iridescent blue and non-iridescent red feathers. Through morphological and chemical analyses, we confirm that hollow, flattened melanosomes in iridescent feathers are eumelanin-based, whereas melanosomes in non-iridescent feathers are solid and amorphous, suggesting that high pheomelanin content underlies red coloration. Intriguingly, the nanoscale arrangement of melanosomes within the barbules was surprisingly similar between feather types. After creating a new genome assembly, we use transcriptomics to show that non-iridescent feather development is associated with genes related to pigmentation, metabolism, and mitochondrial function, suggesting non-iridescent feathers are more energetically expensive to produce than iridescent feathers. However, iridescent feather development is associated with genes related to structural and cellular organization, suggesting that, while nanostructures themselves may passively assemble, barbules and melanosomes may require active organization to give them their shape. Together, our analyses suggest that iridescent feathers form through a combination of passive self-assembly and active processes. 

Hormones mediate physiological and behavioral changes in adults as they transition into reproduction. In this study, we characterize the circulating levels of five key hormones involved in reproduction in rock doves ( Columba livia ): corticosterone, progesterone, estradiol, testosterone, and prolactin using univariate and multivariate approaches. We show similar patterns as previous studies in the overall patterns in circulating levels of these hormones, i.e., testosterone (males) and estradiol (females) high during nest-building or egg-laying, prolactin increasing at mid-incubation and peaking at hatching (both sexes), and elevated corticosterone levels in later incubation and early nestling development. In our investigation of hormone co-variation, we find a strong correlation between prolactin and corticosterone across sampling stages and similarities in earlier (early to mid-incubation) compared to later (late incubation to nestling d9) sampling stages in males and females. Finally, we utilized experimental manipulations to simulate nest loss or altered caregiving lengths to test whether external cues, internal timing, or a combination of these factors contributed most to hormone variation. Following nest loss, we found that both males and females responded to the external cue. Males generally responded quickly following nest loss by increasing circulating testosterone, but this response was muted when nest loss occurred early in reproduction. Similar treatment type, e.g., removal of eggs, clustered similarly in hormone space. These results suggest internal drivers limited male response early in reproduction to nest loss. In contrast, circulating levels of these hormones in females either did not change or decreased following nest manipulation suggesting responsiveness to external drivers, but unlike males, this result suggests that reproductive processes were decreasing. 

Abstract Apoptosis is a fundamental feature of multicellular animals and is best understood in mammals, flies, and nematodes, with the invertebrate models being thought to represent a condition of ancestral simplicity. However, the existence of a leukemia-like cancer in the softshell clam Mya arenaria provides an opportunity to re-evaluate the evolution of the genetic machinery of apoptosis. Here, we report the whole-genome sequence for M. arenaria which we leverage with existing data to test evolutionary hypotheses on the origins of apoptosis in animals. We show that the ancestral bilaterian p53 locus, a master regulator of apoptosis, possessed a complex domain structure, in contrast to that of extant ecdysozoan p53s. Further, ecdysozoan taxa, but not chordates or lophotrochozoans like M. arenaria, show a widespread reduction in apoptosis gene copy number. Finally, phylogenetic exploration of apoptosis gene copy number reveals a striking linkage with p53 domain complexity across species. Our results challenge the current understanding of the evolution of apoptosis and highlight the ancestral complexity of the bilaterian apoptotic tool kit and its subsequent dismantlement during the ecdysozoan radiation. 

Investigation of the negative impacts of stress on reproduction has largely centered around the effects of the adrenal steroid hormone, corticosterone (CORT), and its influence on a system of tissues vital for reproduction—the hypothalamus of the brain, the pituitary gland, and the gonads (the HPG axis). Research on the action of CORT on the HPG axis has predominated the stress and reproductive biology literature, potentially overshadowing other influential mediators. To gain a more complete understanding of how elevated CORT affects transcriptomic activity of the HPG axis, we experimentally examined its role in male and female rock doves ( Columba livia ). We exogenously administrated CORT to mimic circulating levels during the stress response, specifically 30 min of restraint stress, an experimental paradigm known to increase circulating CORT in vertebrates. We examined all changes in transcription within each level of the HPG axis as compared to both restraint-stressed birds and vehicle-injected controls. We also investigated the differential transcriptomic response to CORT and restraint-stress in each sex. We report causal and sex-specific effects of CORT on the HPG transcriptomic stress response. Restraint stress caused 1567 genes to uniquely differentially express while elevated circulating CORT was responsible for the differential expression of 304 genes. Only 108 genes in females and 8 in males differentially expressed in subjects that underwent restraint stress and those who were given exogenous CORT. In response to elevated CORT and restraint-stress, both sexes shared the differential expression of 5 genes, KCNJ5 , CISH , PTGER3 , CEBPD , and ZBTB16 , all located in the pituitary. The known functions of these genes suggest potential influence of elevated CORT on immune function and prolactin synthesis. Gene expression unique to each sex indicated that elevated CORT affected more gene transcription in females than males (78 genes versus 3 genes, respectively). To our knowledge, this is the first study to isolate the role of CORT in HPG genomic transcription during a stress response. We present an extensive and openly accessible view of the role corticosterone in the HPG transcriptomic stress response. Because the HPG system is well conserved across vertebrates, these data have the potential to inspire new therapeutic strategies for reproductive dysregulation in multiple vertebrate systems, including our own.